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Random multifractals occur in particular at critical points of disordered systems. For Anderson 
localization transitions, Mirlin and Evers [PRB 62,7920 (2000)] have proposed the following scenario 
(a) the Inverse Participation Ratios (I.P.R.) Yq{L) display the following fluctuations between the 
disordered samples of linear size L : with respect to the typical value Yq^^{L) = g'"^''^' ~ i^-'^typM 
that involve the typical multifractal spectrum Ttyp{q), the rescaled variable y = Yq(L)/Yq''^{L) is 
distributed with a scale- invariant distribution presenting the power-law tail 1/y^^'^'', so that the 
, disorder-averaged I.P.R. Yq{L) ~ j^-^aviq) jj^yg niultifractal exponents Taviq) that differ from the 

typical ones Ttyp{q) whenever Pq < 1; (b) the tail exponents Pq and the multifractal exponents are 
related by the relation l3qTtyp{q) = Tav{qPq)- Here we show that this scenario can be understood 
by considering the real-space renormalization equations satisfied by the I.P.R. For the simplest 
multifractals described by random cascades, these renormalization equations are formally similar to 
' the recursion relations for disordered models defined on Cayley trees and they admit travelling-wave 

solutions for the variable (IriYq) in the effective time t^ff = InL : the exponent rtyp{q) represents 
the velocity, whereas the tail exponent Pq represents the usual exponential decay of the travelling- 
pH ' wave tail. In addition, we obtain that the relation (b) above can be obtained as a self-consistency 

f— I I condition from the self-similarity of the multifractal spectrum at all scales. Our conclusion is thus 

I , that the Mirlin-Evers scenario should apply to other types of random critical points, and even to 

^ ' random multifractals occurring in other fields. 

' ^ . 

Ti : 

a • I. INTRODUCTION 

^ ^ ^ To explain the motivations of the present work, we first need to recall how the ideas of multifractality on one hand, 

(— ( . and the idea of travelling waves on the other hand, have turned out to play a role in the field of disordered systems. 

o ■ 
o 



o 
o 



Multifractality and critical disordered systems in finite dimension 



^ ' Multifractality is a notion that has first appeared in fiuid dynamics to characterize the statistical properties of 
, turbulence (see the book [l| and references therein). Among the various areas where the multifractal formalism has 
• then turned out to be relevant (see for instance p-Q and references therein) , the case of critical points in the presence 
of frozen disorder is of particular interest. The multifractal character of critical eigenfunctions at quantum Anderson 
. . localization transitions has been the subject of very detailed studies (see the reviews d, and references therein). 
' More generally, multifractality of order parameters and correlation functions is expected to be a generic property 
of random critical points whenever disorder is relevant : it has been found in particular in disordered classical spin 
models like random ferromagnets spin-glasses or random field spin systems |18| - |20H. as well as in disordered 

polymer models like directed polymers in random media [2l| or disordered wetting models [2^ . The only exceptions to 
these multifractal behaviors seem to be the "multiscaling" behaviors [23j , which are even stronger than multifractality, 
that have been found for some critical correlation functions in disordered quantum spin-chains governed by "Infinite 
disorder fixed points" [2^ . 



Travelling waves and disordered systems 



Localized waves that propagate by keeping a fixed shape, have been first discovered in fluid dynamics in 1834 by 
J.S. Russel who wrote : " that singular and beautiful phenomenon which I have called the Wave of Translation" [23] ■ 
The name " Wave of Translation" has not survived, but the idea has flourished under other names. 'Solitary waves' 
or 'solitons' have been found in many areas of physics where non-linear equations of motion occur (see the book [2^ 
and references therein). 'Travelling waves' also appear in particular in the context of front propagation into unstable 
states [i^l : in many cases called "pulled fronts" , the velocity is actually determined by the ex pon entially small tail 
invading the unstable state, and can be thus determined by a linear analysis in the tail region |27l. [28j . It turns out 
that this type of travelling waves also appear in the field of disordered systems, but with of course different meanings 
for the space and time variables with respect to usual spatio-temporal waves. It is useful to distinguish three cases : 
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(i) In disordered models defined on Cayley trees , it is the probability Pi(A) of some observable A that propagates 
without deformation in the effective 'time' corresponding to the length L along the tree 

teff = L (1) 

This property was discovered by Derrida and Spohn (29j on the specific example of the directed polymer in a random 
medium^where the observable A of interest is the free-energy, and was then found in various other statistical physics 
models [30| . This travelling- wave propagation of probability distributions have also been found in quantum models 
defined on Cayley trees, in particular in the Anderson localization problem [3ll - l35j | and in some superfluid-insulator 
transition [36| . The conclusion is thus that the recursion relations that can be written for observables of disordered 
models defined on trees naturally lead to the travelling wave propagation of the corresponding probability distributions. 
This property is not limited to the discrete Cayley trees, but actually still holds for continuously branching trees pot - 

(ii) For some two-dimensional disordered models , travelling waves in the effective time given by the logarithm of 
the system size L 

teff ^ In L (2) 

have been found, first for Dirac fermions in a random magnetic field by Chaman, Mudry and Wen [stI [38|. and then 
for disordered XY models by Carpentier and Le Doussal [39j , who have derived non-linear renormalization equations 
that admit travelling wave solutions. This approach has been then used to study other two-dimensional related models 
[40| , as well as the problem of a particle in a logarithmically-correlated random Gaussian potential in finite dimension. 
For this type of random energy models with logarithmically correlated potentials, many recent developments can be 
found in (4l| - |43| . 

(iii) In finite- dimensional Anderson localization models exactly at criticality , where the eigenfunctions become 
multifractal (see the reviews [l^l or section III Al below) , Evers and Mirlin 44] have proposed that the probability 
distribution PL{lnYq) of the logarithm of the Inverse Participation Ratios Yq (which are the order parameters of the 
Anderson transition) propagate as travelling waves in the effective time tef / = In i given again by the logarithm of 
the system size L. This property has been checked numerically in various Anderson localization models in various 
dimensions [13, [4^ , and has been recently obtained by a functional renormalization method in dimension d = 2 -\- e 
by Foster, Ryu and Ludwig j46j . 

In summary, travelling- wave propagation of probability distributions have been found (i) in most disordered models 
defined on Cayley trees , (ii) in some specific two-dimensional disordered models or for a particle in a logarithmically- 
correlated random potential, and (iii) at Anderson localization transitions in finite dimension. 



C. Multifractality and travelling waves 

We believe that the case (iii) described above, concerning Anderson localization transitions, should actually apply 
to all random critical points in any finite dimension that are characterized by multifractal properties. For the directed 
polymer in a random medium of dimension I -I- 3, we have indeed found numerically the presence of travelling waves 
at criticality [2l|. More generally, besides random critical points, we propose that random multifractal measures are 
generically related to the travelling wave propagation, in the effective time t — InL, of the probability distributions 
of the I.P.R. associated to the multifractal measure. Since the case of arbitrary multifractals is clearly beyond the 
scope of this paper, we will restrict our analysis here to the simplest multifractal measures, namely the random 
cascade models that have been much studied in the context of turbulence [l| . These models satisfy simple real-space 
renormalization equations that are formally similar to the recursion relations for disordered models defined on Cayley 
trees (see case (i) described above). This interpretation thus allows to understand the presence of travelling- wave at 
criticality in finite dimension, as a consequence of a hierarchical real space renormalization procedure on a appropriate 
tree structure in the renormalization scale InL. 

The paper is organized as follows. In section |lTl we recall the multifractal notations and the Evers-Mirlin scenario 
concerning the travelling waves that occur at Anderson localization transitions. In section IIIIl we write real-space 
renormalization equations to describe how the I.P.R. evolve upon coarse-graining. In section IIVI we show that for 
random cascade models, these renormalization equations for the probability distributions of the I.P.R. admit travelling- 
wave solutions. In section |Vl we obtain that the self-similarity of the multifractal spectrum at all scales imposes the 
Mirlin- Evers relation PqTtyp{q) — TaviqPq) that relates the tail exponents (3q of the travelling wave, to the typical and 
disorder-averaged multifractal spectra iTfyp{q), Tav{q))- Our conclusions are summarized in section IVll 
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II. REMINDER ON THE TRAVELLING WAVES AT ANDERSON LOCALIZATION TRANSITIONS 
A. Reminder on typical and averaged multifractal spectra 

At Anderson localization transitions, it is convenient to associate to a critical eigenstate ipLi'r) defined on a volume 
L'^ the measure 

A^L(r) = |VL(r)P (3) 

normalized to 

/ d%i(r) = l (4) 

The Inverse Participation Ratios (I.P.R.) of arbitrary order q are defined by 

Y,{L) ^ I d^vf^lir) (5) 

As a consequence of the normalization of Eq. 01 one has the identity Yq=i{L) — 1. The localization/delocalization 
transition can be characterized by the asymptotic behavior in the limit L — > oo of the Yq{L). In the localized phase, 
these moments Yq{L) converge to finite values 

yJ-(L = oo) > (6) 
In the delocalized phase, the decay of the moments follows the scaling 

r„^'='°'=(i) oc L-(«-i''^ (7) 

Exactly at criticality, the typical decay of the Yq{L) defines a series of generalized exponents Ttyp{q) ~ {q~ l)Dtyp{q) 



where the notation A denotes the average of the observable A over the disordered samples. The notion of multifractality 
corresponds to the case where Dtyp{q) depends on q, whereas monofractality corresponds to Dtyp{q) = est as in Eq. 
([7]). The exponents Dtyp{q) represent generalized dimensions [H : Dtyp{0) represent the dimension of the support of 
the measure, here it is simply given by the space dimension Dtyp{0) = d; Dtyp{l) is usually called the information 
dimension Q , because it describes the behavior of the 'information' entropy 

SL = -y^UL{r}^nnL{r) ^ ~dqYq{L)\q=i cx Dtyp{l)\nL (9) 

r 

Finally Dtyp{2) is called the correlation dimension 2] and describes the decay of 



L— ^cxj 

In the multifractal formalism, one also introduces the singularity spectrum ftypiot) defined as follows : in a sample 
of size L'^, the number Mhia) of points r where the weight scales as L~" behaves typically as 

M^i^^ia) oc L-^*«^(") (11) 

L— >oo 

The saddle-point calculus in a of the I.P.R. 

Y^yPiL) ~ y" da L/'«p(") (12) 

yields the usual Legendre transform formula 

- Ttyp{q) = max [Jtyp{a) - qa\ (13) 
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Following 0, many authors consider that the singularity spectrum has a meaning only for ftyp{a) > [47l - [50| . 
However, when multifractality arises in random systems, disorder-averaged values may involve other generalized 
exponents [sil - ls^ than the typical values (see Eq. |H]), and it is thus useful to introduce another series of generalized 
exponents Tav{q) = (g - l)Dav{q) [S] 

cx = L~(«-i)^""(«) (14) 

For these disorder averaged values, the corresponding singularity spectrum favict) may become negative fav{a) < 
[3. [5ll - l55| to describe rare events (cf Eq. [TT]) . The difference between the two generalized exponents sets Dtypio) and 
DavW) associated to typical and averaged values has for origin the broad distribution of the I.P.R. over the samples 
1441 1451 as we now describe. 



B. Probability distributions of the I.P.R. Yq{L) over the samples 

The scenario proposed by Evers and Mirlin [3, EH] in the context of quantum localization models is as follows : the 
probability distribution of the logarithm of the Inverse Participation Ratios of Eq. [5] becomes scale invariant around 
its typical value 4445|], i.e. 



lnr,(i) = lny,(L) + = \nYlyP{L) + Uq (15) 

where Uq remains a random variable of order 0(1) in the limit L — > oo. According to [43 | the probability distribution 
Ghi'^q) generically develops an exponential tail 

Gooiuq) oc e-P-^- (16) 



As a consequence, the ratio yq ~ Yq{L) / Y^^p {L) — e"' with respect to the typical value YqyP{L) = e'"^'*^-^) presents 
the power-law decay 

The conclusions of [3 then as follows : 

(i) The typical singularity spectrum has a meaning only for ftyp{a) > 0, i.e. it usually exists only on a finite interval 

where the termination points a± satisfy ftyp{a.-) = = ftyp{a+). Denoting q± the corresponding values of 
q by the Legendre transformation of Eq. I13[ one obtains the linear behaviors outside the interval q^] 

Ttypiq) = qa- for q < q- (18) 
Ttypiq) = qa+ for q > q+ (19) 

(ii) The disorder-averaged singularity spectrum fav{ct) has a meaning outside this interval where favict) < and 
describes the probabilities of rare events. In the region where both exist, they are expected to coincide 

fav{a) = ftyp{a) if a+ < a < a- (20) 

Equivalently, the disorder-averaged exponents Tav{q) are expected to coincide with the typical exponents Ttyp{q) on 
the interval 

Tav{q) = nyp{q) if q- <q<q+ (21) 

whereas Tav{q) will be different from Eq. [12] outside the interval 

(iii) from the point of view of the tail exponents (3q of Eq. [T71 this means that 

pq > 1 if (7_ < (? < g+ (22) 

Pq <l li q < q- ov q> q+ (23) 

/3g± = 1 (24) 

(iv) finally, Mirlin and Evers [i^ have derived in special cases the following relation 

PqTtypiq) ^ TaviqPq) (25) 

and they have conjectured its generic vahdity [13, Hi- 

This scenario has been tested numerically for various Anderson transitions (see the review [To| and the more recent 
works [55*1). Wc have found previously that this scenario also applies to the transition of the directed polymer in 
dimension 1 + 3 [21.] . In the following, we justify this scenario via a real-space renormalization analysis. 
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III. REAL-SPACE RENORMALIZATION ANALYSIS 



In critical phenomena, it is well known that critical properties are stable under coarse-graining. This explains their 
universal character (independence with respect to microscopic details) and why renormalization is an appropriate 
framework. Similarly for random critical points, the multifractal spectrum is expected to be stable under coarse- 
graining Q. It is thus natural to consider how the Inverse Participation Ratios (I.P.R.) evolve upon coarse-graining. 



A. Evolution of the I.P.R. upon coarse-graining 



To go from the microscopic scale I — 1 to the macroscopic scale I — L oi the whole system, it is convenient to 
introduce intermediate scales Im regularly placed on a logarithmic scale as follows. For definiteness, we consider a 
discrete system with L = b^^ and introduce the intermediate scales 

l^ = ir with m = 0,l,..,M (26) 

Then the whole volume L'^ = ij^j can be decomposed into 

K = h'^ (27) 

subvolumes of sizes Im_i, and the process can be recursively iterated : each volume of size 1"^+! = Klf^ is made of K 
volumes of sizes ZJ^, denoted here by an index j = 1,2, ..K. We consider a (non-normalized) positive field ^Ji{r)^ and 
associate to each volume (i) of size 1^. the integrals of ^.''{r) over the corresponding volume 

Z«(m) = j d'^rn'^ir) (28) 
The corresponding I.P.R. are then defined by the ratios 



^zW(m) 

Upon coarse-graining, the integrals Zq are simply additive for any q 

K 

Z,(m + l) = ^ZW(,n) (30) 

i=l 

whereas the I.P.R. satisfy the recursions 

K 

Yq{m + l) = Y. Y^^{m) (31) 

1=1 

where the coefficients are given by the weights 

-.(m) - (32) 

that represent the ratios of the normalisation in the volume (i) of size with respect to the normalisation of the 
volume of size Im+i- By construction one has the following constraint 

K 

Y,w,{m)^\ (33) 

1=1 

For instance for Anderson localization models where the additive positive field /i(r) is given by Eq. [31 one expects 
that at sufficiently large scale 

(i) in the delocalized phase, the K weights Wi all converge towards the same value 1/ K, so that the system becomes 
asymptotically homogeneous at sufficiently large scales. 

(ii) in the localized phase, one single weight converge to 1, whereas all other {K — 1) weights converge to zero. 

(iii) exactly at criticality, the K weights remain finite in contrast to (ii), but they remain distributed with a non 
trivial distribution, in contrast to (i). 
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B. Notion of random cascade models 



Since the general case where the weights are correlated among generations is more difficult to analyse, we will 
consider from now on the much simpler case where the weights of different generations are uncorrelated, and where 
the non- negative weights {wi, wk) corresponding to an elementary coarse- graining step in Eq. |3T]are drawn with 
some fixed probability distribution QKwi, ...,wk), independent of the generation m, symmetric in its K arguments, 
and satisfying only the normalization constraint of Eq. [33] : 

QI{wi,...,wk) = RI{wi,...,wk)S Wi(m) - 1^ (34) 

This type of random cascade models has been much studied in the context of turbulence to describe the spatial 
distribution of energy dissipation (see the book [l| and references therein). In this context, the weights {wi) are 
usually called 'cascade generators' or 'multipliers' or 'breakdown coefficients'. Various forms have been proposed over 
the years, in particular log- normal [iH, bimodal [stJ, log-stable [H, HI], log-Poisson [60l - [6^ . log- infinitely-divisible 



[63| : see [6j| for a comparative test of these various cascade generators. 

However besides the specific form of the distribution of these weights, the important hypothesis is of course the 
independence of the weights of different generations. For turbulence, the validity of this hypothesis is discussed 
in [gsI |66| . For random critical points, this hypothesis is not expected to be valid, but one expects instead some 



Markovian structure 67[. Nevertheless, since random cascade models are clearly the simplest multifractal models, 
it is important to understand in detail their properties. In the following, we show that travelling-waves appear very 
naturally in random cascade models, and that their properties are agree with the Mirlin-Evers scenario concerning 
Anderson transitions. 



IV. TRAVELLING- WAVE ANALYSIS FOR RANDOM CASCADE MODELS 



In the present section, we analyse the real-space renormalization Eq. [31] for the case of random cascade models 
characterized by some fixed distribution Ql{wi, ...,wk) (see Eq. I34|) . The renormalization Eq. [3l]is then analogous 
to recursion equations for disordered models defined on a Cayley tree of branching number K, and it is thus natural 
to obtain travelling-wave propagation of probability distribution as for the directed polymer model [1^ mentioned in 
the introduction. The fact that renormalization of multifractals in finite dimension involves a hierarchical structure 
analogous to the Cayley tree has been already stressed in various contexts [37,-.40> .68] . 



A. Travelling wave Ansatz for the Inverse Participation Ratios (I.P.R.) 



We look for solutions of the real-space renormalization Eq. [21] with the following Ansatz 

Y,im) = e--^"'y, (35) 

where Vq is a constant depending on q, and where yq is a random variable whose distribution does not depend upon 
the generation m. In logarithmic scale, this corresponds to the travelling wave form 

hiYq{m) — —Vqm + Iny, (36) 

where Vq represents the velocity with respect to the variable m = (In Z^)/ (In 6). The velocity Vq is thus directly related 
to the typical exponent Ttyp (q) of Eq. M 

Vq = {\nb)Ttyp{q) (37) 
With the Ansatz of Eq. [HS] the renormalization Eq. [21] becomes 

K 

yq = e"' t"'^]' y^^^ (38) 
1=1 

for the reduced random variables yq. The probability distribution Pq{yq) should be stable by this iteration 

P;{yq) = I dw,...dwKQl{wu ....wk) I dy«P;(y«)...y' dyf ^^^(yf ))<5 U - e"' (^9) 
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B. Tail analysis 

A gen eric property of multiplicative stochastic processes is to lead to probability distribution presenting power-law 
tails |69l - [74| . Here for the specific case of random multiplicative cascade models, one expects also that the stable 
distribution Pqiuq) solution of Eq. 131] will present a power-law tail in j/^ 

(40) 

where the exponent I3q is not fixed for the moment, and satisfies only the condition /3g > to have a normalizable 
probability distribution. In the travelling wave language of Eq. 1361 this power-law tail is equivalent to the exponential 
tail e"^'"' for the variable u„ = In Yq(m) -I- mvq —Inyq. So this corresponds exactly to the standard exponential tail 
analysis of travelling fronts [13, • 

The stability of the power- law tail of EqHHliii the region yq — > -|-cxd means that at leading order, only one of the K 
variable y^'^ in Eq. [51] becomes large with a probability also given by the tail of Eq. [30] : after the introduction of a 
factor K to choose one of the K variables j/g*-* , we may assume the choice i = 1 leading to 



yq J J J J (ijq ) Pq 

K [ dwi...dwKQl{wi,...,WK)wf'e''^^^ (41) 



yq -J 

One thus obtain the following compatibility equation 



1 = Ke^'^^^wf" (42) 
in terms of the partial moment of the joint probability distribution Ql{wi, ...,wk) of Eq. 1341 

dwidw2--dwKQb{wi, ...,wk)wI (43) 

Eq. [42] means that each mode /3 is characterized by the velocity Vq (/3) given by 

Vq{P)^~-^\n(K^) (44) 



C. Selection of the tail exponent jiq and of the velocity Vq of the travelling wave 

In the field of travelling waves, the selection of the tail exponent /? of Eq. [30] and of the corresponding velocity 
v{j3) of Eq [36] usually depend on the form of the initial condition |27l - l29| . In our present case, the initial condition 
is completely localized YJrn = 0) = 1 at the lowest scale lm=o = 1- One then expects that the solution that will be 
dynamically selected |27l - [29| corresponds to the tail exponent /S^'^^'^'^ and to the velocity v^'^^'^'^ = Vq{l3q'^^'^'^) determined 
by the following extremization criterion 



(45) 

f^selec 

In summary, for any given distribution (5^(wi, wj<-) that defines a random cascade model (see Eq. I34[) . the 
properties of the travelling waves can be obtained as follows : one has to compute the partial moment of Eq. 1431 and 
to solve Eq. [45] in order to determine the tail exponent jiq and the velocity Vq that are dynamically selected. The 
typical multifractal exponents Ttyp{q) are then obtained from the selected velocities by Eq. [37] 



= [d,vqmp=p..... = In (Kwf) - In 



ySelec 

^'y^^'^^ = (hTfo) 



(46) 
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V. RELATIONS BETWEEN THE TAIL EXPONENTS AND THE MULTIFRACTAL EXPONENTS 

As explained in the previous section, the selection criterion of Eq. |45] is usually the final outcome of a travelling 
wave analysis. However here in our real-space renormalization framework, we still have some freedom in the choice 
of the rescaling factor b introduced in Eq. [521 Of course, the final multifractal exponents Ttyp{q),Taviq) should not 
depend on the choice of the coarse-graining scale b. In the present section, we use this freedom to consider the case 
of large b, and we show that some self-consistency conditions arise. 



A. Self-consistency conditions when the rescaling factor b is large 

When the rescaling factor b introduced in Eq. [55] is large, the number K — b"^ oi buildings blocks (see Eq. \T7]\ 
of a single renormalization step also becomes large. Then from the self-similarity of the multifractal spectrum at all 
scales [91], one obtains that the elementary weights Wi should themselves follow the multifractal statistics, i.e. their 
associated I.P.R. 

y,{b) ^Y.""' (47) 



i=l 



should have for typical scalings 



and for averaged scalings 



yiy^ij}) cx fe-^'^ptrt (48) 

^ 6— >4-oo 



yp{b)^Kw\ oc fe---(p) (49) 

in terms of the typical and averaged multifractals exponents Ttyp{p) and Tav{p) introduced in Eqs l8l and 1141 The 
velocity of Eq. |44] then reads at leading order for large b 

v,{(3) = -iln = l2^fcMin6 (50) 

Taking into account Eq. I46[ we thus obtain the following consistency relation 

selec I oselec\ 

that relates the typical exponent Ttyp{q) to the selected tail exponent /J^'^''^'^ and to the averaged multifractal exponent 
Tav{p = qP^^ff^)- As recalled around Eq. [25l the relation of Eq. [51] has already been derived in special cases by Mirlin 
and Evers (44| and has been conjectured to be general [IB]. Our present derivation from a self-consistency condition 
of the travelling wave analysis is thus in favor of the general validity of this formula. 

The relations of Eq. [51] introduces non-trivial constraints on the multifractal exponents, which are not always 
compatible with the usual selection criterion of Eq. [45] based on the extremization of the velocity. Let us first describe 
an explicit case before returning to the general case. 



B. Example of Gaussian multifractality 

The simplest multifractal spectrum corresponds to the following Gaussian form for the disorder-averaged multifractal 
spectrum 



^{q) = d{q-l)\l-^^j (52) 

In particular, this Gaussian form s ap pears in various models in the weak multifractality regime [lo| . in particular in 
perturbation theory in d = 2 -|- e [46[ . 
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1. Consequences of the self-consistency Eq. \51\ 
In the region where T(yp"'*'*(g) — T^y^^'^{q) corresponding to > 1, Equation [5T] yields [lOj 

P. = f. (53) 
This solution is consistent for /3g > 1, i.e. in the interval 

-qc<q<+qc (54) 

In the region where r^p'^'^^iq) ^ t^v"'^'''' il) corresponding to I3q < 1, the typical exponents vary linearly with q 

r?yr'{q) ^qd{l- for |g| > (55) 



Equation [5T] then yields \\^\^ 

q 



= n (56) 



2. Analysis via the selection criterion of Eq, \51\ 
From Eq. [SH we obtain the velocity as a function of the tail exponent /3 using Eq. (SD] 

"™(^-^ln6 = ^-^diPq - 1) fl - ^'^ 

P fi \ q 



v^iP) = ^i^^A^i^lnfo = S^di^q - 1) ( 1 - ^ ) (57) 



Its derivative 

df,vg{(3) ^ d{\nb) 



' 1 g2 

yields the following solution (3'^^^^'^ > for the selection criterion of Eq, 

elec 9' 



(58) 



and the corresponding selected velocity reads 



V 



,^selec 
9 



^ _ sgn(g) \ 



z;,(/3-'-)-(ln6)dg 1-^^ (60) 



qc / 



The typical multifractal exponent reads (Eq. |46 



/ sgn(g)y 

i.e. one recovers the correct result only in the linear regions of the typical spectrum, where it is different from the 
disorder-averaged spectrum . This seems to indicate that in the region where Ttyp{q) = Tav{q), this additional constraint 
'''typiq) = Taviq) fixes completely the tail exponent and the velocity to values that do not satisfy the usual selection 
criterion of Eq. 1511 This phenomenon seems generic even beyond the Gaussian case as we now explain. 

C. General case : competition with the usual velocity selection 

For the general case, the selection condition of Eq. |35] for (3^^''^'' becomes 

O=[qPTUqP)-rav{qP)]0s.u. (62) 
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In terms of the solution p* of the foUowing equation 



= K.(p)-^a.(p)],=p. (63) 



(64) 



the selected tail exponent reads 

and the corresponding velocity is given by 

vf^" = vJpf'"') = °"^^ 7 In 6 (65) 

y y v y / j^selec ^ ' 



Eq. [571 yields the following typical exponent 

ln& " |p*|sgn(9) 



which is linear in q. 

In summary, as in the Gaussian case, we obtain that 

(i) the usual selection criterion of Eq. |45] is compatible with the self-consistency Eq. [51] only in the region where 
the typical and averaged spectra differ Ttyp{q) ^ Taviq) where /Sq < 1. 

(ii) In the region where Ttyp{q) — Tav{q) and (3q > 1, the self-consistency Eq. [51] indicates that the selected tail 
exponent is the solution /3^^'^'^ > 1 of 

f^f'^ra.iq) = Ta.iqPfn (67) 

(the other trivial solution of this equation being (3 — 1), and does not correspond to the usual selection criterion of 
Eq. M 

VI. CONCLUSIONS AND PERSPECTIVES 

In this paper, we have proposed that the Mirlin-Evers scenario concerning the probability distributions of In- 
verse Participation Ratios (I.P.R.) at Anderson transitions can be better understood by considering the real-space 
renormalization equations satisfied by I.P.R. upon coarse-graining. For the simplest multifractal models, namely the 
random cascade models, we have shown that these renormalization equations are formally similar to the recursions 
equations for disordered models defined on Cayley trees, and that, as a consequence, they admit travelling wave 
solutions of the type known as "pulled fronts" , where the velocity is actually determined by a linear analysis in the 
exponentially-small tail region. Finally, we have obtained that the self-similarity of the multifractal spectrum at all 
scales imposes the Mirlin-Evers relation of Eq. [2S] that relates the tail exponents /3q of the travelling waves, to the 
typical and disorder-averaged multifractal spectra {Ttyp{q),Tav{q))- This shows that random cascade models already 
capture many properties that have been previously found at Anderson transitions. It would be thus interesting in 
the future to show that the travelling wave solutions of the renormalization equations persist beyond random cascade 
models, since at random critical points, one expects that the strict statistical independence of the weights of different 
generations is not valid, but one expects instead some Markovian structure (67j . Further work is needed to formulate 
correctly appropriate models of Markovian cascades. 

Nevertheless, since random cascade models are clearly the simplest multifractal models, and since their properties 
are very similar to the properties found previously for the more complex case of Anderson transitions, we believe 
that our present real-space renormalization analysis is in favor of a wide validity of the Mirlin-Evers scenario beyond 
Anderson transitions : 

(i) the first generalization concerns all random critical points in finite dimension, where the local order parameter and 
the correlation functions generically display multifractal statistics (see the discussion in section [T^ of the Introduction) . 
For instance, we have found numerically the presence of travelling waves at criticality for the directed polymer in a 



random medium of dimension 1-1-3, [21|. For other many-body random phase transitions like classical disordered 
spin models, these travelling wave properties should also appear after an appropriate translation (see [13] for such a 
translation between multifractality at Anderson transitions and in classical disordered spin models) . It seems that for 
two-dimensional disordered Potts models where multifractality has been much studied [TT| - [l7l [67j only the disorder- 
averaged multifractal spectrum has been considered up to now. It would thus be very interesting to study numerically 
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the probability distribution over the disordered samples of the appropriate observables to test the travelhng-wave 
scenario and to measure the corresponding tail exponents Pg. Some probability distributions over the disordered 
sainples of extensive observables like the susceptibility have been already measured for some random critical points 
in [Ta - lTSj , but a quantitative analysis of the tails remains to be done. 




(ii) besides phase transitions in disordered systems, a further generalization concerns other areas of physics where 
multifractality occurs. We believe that the present real-space renormalization analysis on the multifractal measure 
keeps its validity, provided the notion of fluctuations between different realizations of the multifractal cascade has a 
physical meaning. Indeed, in disordered systems, each disordered sample is characterized by a given realization of 
the multifractal cascade, so that the fluctuations between cascade realizations correspond to the sample-to-sample 
fluctuations, which play in major role in the understanding of disordered systems. In other fields where multifractality 
occurs, one should first clarify the physical meaning of a given realization of the cascade to see whether it is interesting 
to distinguish between typical exponents and disorder-averaged exponents, and to introduce probability distributions 
over cascade realizations. 



It is a pleasure to thank Olivier Giraud for a very useful discussion on multifractality and tail exponents. 
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